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We present a numerical and analytical study of scalar <fr theory at finite temperature 

. with a renormalized 2-loop truncation of the 2PI effective action. 
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1. Introduction 

^ : 

Q-i! The thermodynamical properties of a physical system can be extracted from the effective 
Q_i! potential, which is given as a function ~f[<f>] of a condensate field 0. In the thermodynamic 
limit, its stationary point defines the free energy T as 

b3: -F=lim 7 [0], with ^1 =0. (1) 

d; 

Once the free energy is known, all other thermodynamic quantities, such as the entropy 
or energy density, can be easily derived. At finite temperature collective phenomena are 
known to modify substantially the properties of the elementary excitations, so the use 
of a perturbative expansion around the free theory to calculate the effective potential 
becomes questionable and often leads to inconsistencies. One needs then to consider 
nonperturbative schemes. For situations where the effects of the collective phenomena can 
be conveniently captured by a suitable modification of the 2-point functions, resummation 
schemes based on the 2PI effective action can be very useful [|2]. The 2PI effective action 
provides a reorganization of the perturbative expansion around dressed 2-point functions, 
which are determined self-consistently for a given approximation/truncation. In this work 
we study the thermodynamics of scalar </> 4 theory using a 2-loop truncation of the 2PI 
effective action. 

2. 2-loop truncation of the 2PI effective action 

For scalar 4 theory the 2PI effective action can be parametrized as [Q] 

r 2 pi[0, G}= l -<P- G 1 ■ + ^Tr [luG- 1 + (G„ 1 - G' 1 ) ■ G] + T int [0, G], (2) 
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with and G generic 1- and 2-point functions and A ■ B a shorthand notation for the 
convolution of A and B. The term r int contains the interactions and can be written as an 
expansion in terms of two-particle-irreducible (2PI) diagrams. Up to 2-loops it is given 
by 1 



-rw0,G] = -x + 7>c+o X+^^e4^ (3) 




Within the approximation, "physical" 1- and 2-point functions and G are determined 
self-consistenly as the stationary points of the 2PI effective action, i.e. 



5T 



2PI 



5G 



and 1 ' 



G[4>U 5(f) 



_ _ = 0. (4) 

G[4>U 



The stationary conditions (HJ) turn into a set of coupled implicit equations for G and 
which have to be solved in order to calculate any physical quantity. In particular, the 
knowledge of the dressed 2-point function G allows one to calculate the effective potential 
as 7 [0] =T\Z- 1 r 2PI [0,G[0]]. 

One of the main complications that arises when dealing with truncations of the 2PI 
effective action is the fact that 2- and higher n-point functions are not uniquely defined [Ej. 
In particular, a given truncation defines two possible 2-point functions and three 4-point 
functions. One of the 2-point functions is given by the stationary value G, while the other 
is related to the inverse curvature of the effective potential as G~ 1 = TV^ 1 [S 2 j[4>]/5(j) 2 ]. 
The three 4-point functions and their corresponding vertices are given in terms of coupled 
Bethe-Salpeter-like equations [ 4 . Concerning renormalization, the ambiguity in the 
definition of the vertex functions implies that there are more than one counterterm of a 
given type [Ej- A given counterterm is determined by applying a renormalization condition 
to the corresponding vertex. For consistency, identical renormalization conditions are 
applied to all counterterms of the same type. It is simpler to consider renormalization 
conditions applied at a reference temperature T* for which the physical field configuration 
is 0* = 0. This requires that > T c if a critical temperature T c exists. Defining the 
self-energy from Dyson's equation S = G^ 1 — Gq 1 , the renormalized equations for S and 
are 

£(P) = 5ml + ^±^0 2 + f T G{K) - ^0 2 ©(P) (5) 



2 2 J K y ' 2 

"T \ 2 i-T 



(m 2 + 6ml) = ^103 + ±±lh^ ^ G{K) - ^0 jT G(K)Q(K), (6) 
with $x a short-hand notation for the standard sum-integral over momentum K at tem- 



perature T, and O(P) = G(K)G(K + P). The bar in has been omitted in the first 
equation to stress the fact that it can be solved for any value. This is needed, for instance, 
to calculate the effective potential j[4>]. The explicit expressions for the counterterms Sm^, 
5m 2 ,, 5\ , 5X 2 and 6X4, can be found in Ref. [Sj. With those counterterms it can be shown 
that, for any value of T and 0, the results for G, and 7[0] are UV-finite. 



1 The Feynman rules are given by: 



= -A, = G , 



Renormalized Finite Temperature <f) 4 theory from the 2PI Effective Action 



3 




po/m 

Figure 1. Spectral function p(po,0) for A = 1, <p/m — 1 and T^jm = 1 and several 
temperatures. 

3. Numerical Analysis 

The gap and field equations (J^J) and (0) are solved numerically in Minkowski space. 
We split the self-energy into a local and a non-local part as S(P) = + S n '(P). The 
momentum dependent part is expanded in both p and | p | using N Chebyshev polynomials 
as S ni (p , |p|) = c ijTi(po)Tj(\p\). The solution to the gap and field equations is 

obtained by solving the matrix equation for the N 2 coefficients Cjj plus the local terms 
£' and/or 2 by means of a multi-dimensional Newton- Raphson method. This numerical 
algorithm is fairly stable and converges after few iterations. For a given N, the numerical 
solution oscillates around the "full" solution, which is convenient for the calculation of 
integrated quantities, such as the effective potential. Unlike lattice-based methods, this 
technique also allows the use of large (small) UV (IR) cut-offs, which is useful for the study 
of renormalization and/or critical phenomena. Although numerically more expensive, the 
advantage of solving the equations in Minkowski space is that one can obtain directly 
spectral properties of the system without the need of analytic continuation. In particular, 
the spectral function p(po, |p|) can be constructed from the knowledge of the real and 
imaginary parts of the retarded self-energy, which come directly from solving equation 
(j3J) (see Fig. 01 for an example result at several temperatures). The width, the position 
of the quasiparticle pole and the effect of multiparticle contributions can be adequately 
studied with the presented algorithm. 

We can also look at the critical behaviour in a system with broken symmetry. A simple 
quantity to compute is the critical temperature which can be read off the 2-point functions 
at vanishing effective mass M. As we discussed previously there are two possible 2-point 
functions (G and G), for which the corresponding effective masses are (in the symmetric 



phase) 



M 2 (T,A) 



m 2 + S'(T,A), 



(7) 



M(T, A) 



G-\p = 0,T,X) 



50 2 



(p = 0,T,A). 



(8) 



Starting from T* and reducing the temperature, the critical values T c and T c are found 
when M 2 (T C , A) and M(T C , A) vanish (see Fig. El left). We can also vary the UV-cutoff 
in the calculation of the critical temperatures and therefore check the validity of the 
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Figure 2. Left: Critical temperatures for several couplings and T*/m 
cut-off dependence for A = 3. 



5. Right: UV 





Figure 3. Left: field expectation value 2 as a function of T for A = 3 and T+/m 
Right: Effective potential. 
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renormalization procedure (see Fig. El right). We find that both quadratic and logarithmic 
divergences (present respectively in T c and T c ) are properly renormalized and hence the 
renormalization is satisfactory. 

Solving in addition the field equation (0) allows us to study the phase transition (see 
Fig. left). We see the transition is clearly of second order. A similar conclusion is 
reached by computing directly the effective potential j[(f>\ (see Fig. left). We checked 
that 7[0], from which a renormalized free energy and pressure can be extracted, is cutoff 
independent. 
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